Numerical Computations of Self - Similarblow - up Solutions of Thegeneralized
نویسندگان
چکیده
The structure of the blow-up in nite time of a solution of the Generalized Korteweg-de Vries equation arising from a perturbed unstable solitary wave is studied numerically. The computed solution is observed to blow-up in the L 1-norm in nite time by forming a spike of innnite height at x = x and at t = t. Scaled coordinates are introduced to examine the detailed structure of the solution in the immediate neighborhood of the blow-up. The appropriately rescaled solution is observed to converge in these coordinates as t ! t ? , indicating self-similar behavior. A best-t solution w() of the nonlinear ODE satissed by self-similar prooles is computed for the statistical data compiled from this convergence. The asymptotics at 1 of this solution of the ODE are studied, and found to coincide with those of solutions w () of the linearized ODE as ! 1. The self-similar part of the solution is also matched (numerically) to the part of the solution more removed from the blow-up point, showing how rapidly decaying initial data can give rise to self-similar blow-up. Heuristic explanations of how nonlinearity and dispersion cooperate to yield existence of a solution w() of the ODE with the desired asymptotics as ! 1 are discussed.
منابع مشابه
NUMERICAL SOLUTIONS OF SECOND ORDER BOUNDARY VALUE PROBLEM BY USING HYPERBOLIC UNIFORM B-SPLINES OF ORDER 4
In this paper, using the hyperbolic uniform spline of order 4 we develop the classes of methods for the numerical solution of second order boundary value problems (2VBP) with Dirichlet, Neumann and Cauchy types boundary conditions. The second derivativeis approximated by the three-point central difference scheme. The approximate results, obtained by the proposed method, confirm theconvergence o...
متن کاملNON-POLYNOMIAL SPLINE SOLUTIONS FOR SPECIAL NONLINEAR FOURTH-ORDER BOUNDARY VALUE PROBLEMS
We present a sixth-order non-polynomial spline method for the solutions of two-point nonlinear boundary value problem u(4)+f(x,u)=0, u(a)=α1, u''(a)= α2, u(b)= β1,u''(b)= β2, in off step points. Numerical method of sixth-order with end conditions of the order 6 is derived. The convergence analysis of the method has been discussed. Numerical examples are presented to illustrate the applications ...
متن کاملPrecise computations of chemotactic collapse using moving mesh methods
We consider the problem of computing blow-up solutions of chemotaxis systems, or the so-called chemotactic collapse. In two spatial dimensions, such solutions can have approximate self-similar behaviour, which can be very challenging to verify in numerical simulations [cf. Betterton and Brenner, Collapsing bacterial cylinders, Phys. Rev. E 64 (2001) 061904]. We analyse a dynamic (scale-invarian...
متن کاملThe comparison of two high-order semi-discrete central schemes for solving hyperbolic conservation laws
This work presents two high-order, semi-discrete, central-upwind schemes for computing approximate solutions of 1D systems of conservation laws. We propose a central weighted essentially non-oscillatory (CWENO) reconstruction, also we apply a fourth-order reconstruction proposed by Peer et al., and afterwards, we combine these reconstructions with a semi-discrete central-upwind numerical flux ...
متن کاملTransition to Blow-up in a Reaction-Diffusion Model with Localized Spike Solutions
For certain singularly perturbed two-component reaction-diffusion (RD) systems, the bifurcation diagram of steady-state spike solutions is characterized by a saddle-node behavior in terms of some parameter β in the system. For some such systems, such as the Gray-Scott model, a spike self-replication behavior is observed as a parameter varies across the saddle-node point. We demonstrate and anal...
متن کامل